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Abstract. I review recent progress in using twistor-inspired methods to compute perturbative 
scattering amplitudes in gauge theory, for application to collider physics. 


INTRODUCTION 

In two years, a new window will open into physios at the shortest distance scales. The 
Large Hadron Collider (LHC) will begin operation at CERN, providing proton-proton 
collisions at 14 TeV center-of-mass energy, seven times greater than the 2 TeV currently 
available inpp collisions at Fermilab’s Tevatron. The LHC luminosity should be a factor 
of 10 to 100 greater than the Tevatron’s. The combined rise in energy and luminosity 
will lead to a huge increase in the production of particles with masses in the range 100- 
1000 GeV, including electroweak vector bosons, top quarks, Higgs bosons, and of course 
new particles, representing physics beyond the Standard Model. 

There are a lot of ideas for physics beyond the Standard Model, many associated with 
the puzzle of electroweak symmetry breaking, and with resolutions of the hierarchy 
problem — why the weak scale is so much smaller than the Planck scale. Supersym¬ 
metry, for example, predicts a host of new particles in the 100-1000 GeV mass range, 
including (in most versions) a stable dark matter candidate. However, many other the¬ 
ories — new dimensions of space-time, new forces, etc. — often make qualitatively 
similar predictions. How can we sort out the predictions of these theories from each 
other, and from the omnipresent Standard Model background at a hadron collider? 

The short answer is that a thorough, quantitative understanding of both the new 
physics signals and the Standard Model backgrounds is required. Much work has gone 
into these problems, stretching back over many decades. This talk will focus on some 
recent developments, novel methods to help compute the backgrounds in particular, that 
have emerged since the Fall of 2003, when Witten introduced twistor string theory and 
explained its relevance to perturbative QCD [Ql]. In truth, the new methods have not 
yet had a direct phenomenological impact, in terms of producing more accurate cross 
sections that have not previously been obtained in any other way. But they have a lot of 
promise, and it should not be long before they do so. 

What are some generic properties of the new physics signals? Except for stable. 


^ Talk presented at the International Europhysics Conference on High Energy Physics, Lisbon, Portugal, 
July, 2005. Research supported by the US Department of Energy under contract DE-AC02-76SE00515. 





FIGURE 1. Typical cascade decay of a gluino to two quarks, a gluon, a W boson which decays to two 
leptons, and a neutralino. Aside from the neutralino, all the final-state particles are essentially massless. 


neutral dark matter candidates, the new massive particles typically decay into “old” 
Standard Model particles: quarks, gluons, charged leptons and neutrinos, photons, Vks 
and Zs. For example, in supersymmetry the superpartner of the gluon, the gluino, may 
be among the heavier superparticles, yet still be copiously produced at the LHC, due 
to its large adjoint color charge. Figure U shows a typical decay cascade, initiated by 
one of the two gluinos (g) in a pair-production event. The quarks and gluons emerge 
as jets of hadrons. The lightest superpartner, a neutralino (X^), is stable and escapes 
the detector. The kinematic signatures of such events are not always clean: There can 
be a large number of observed particles (charged leptons or jets), and no invariant- 
mass bumps, because of the escaping neutralinos and possibly neutrinos (although there 
can be kinematic edges). The escaping neutralinos provide a missing transverse energy 
signal, but Standard Model production of Z bosons, followed by Z decays to neutrinos, 
can mimic this to some degree. 

In order to maximize the potential for the discovery and interpretation of new physics 
at the LHC, we need to quantify the Standard Model backgrounds for processes that 
may contain several jets and (perhaps) a few electroweak bosons. These processes are 
complex, so we should try to take into account any simplifying features. Notice that 
the masses of the observed final-state particles in these reactions (e.g. in fig. [T]) are 
generally negligibly small in these reactions, except for the cases of the W, Z, or top 
quark. Even these particles immediately decay to essentially massless quarks or leptons, 
however. So if we include the decay processes in the description of the event, every 
final-state particle is approximately massless. We can also (usually) neglect the masses 
of the colliding partons (quarks and gluons). In general, then, the backgrounds (and 
many signals) require a detailed understanding of scattering amplitudes for many ultra- 
relativistic (massless) particles - especially the quarks and gluons of QCD. 

Asymptotic freedom ||2] allows us to compute such scattering amplitudes as a pertur¬ 
bative expansion in the strong coupling constant aj(/i), evaluated at a large momentum 
scale /i where it is small. For typical collider processes, /i could be of order 100-200 
GeV, for which (Xs{}l) ~ 0.1. One might expect that the leading-order terms in the ex¬ 
pansion (tree amplitudes) would suffice to get a 10% uncertainty. However, this is not 
the case for hadron collider cross sections; typical corrections from the next-to-leading 
order (NLO) terms in the Us expansion are 30% to 100%. There are several possible rea¬ 
sons for the large corrections, depending on the process: there may be different scales 
involved, leading to large logarithms of the ratio(s) of scales; new partonic subprocesses 
may first arise at NLO; the lowest-order process may have several factors of Us in it; 




and so on. In any event, a quantitative description of collider events requires evaluation 
of cross sections at NLO in QCD, which in turn requires, as input, one-loop amplitudes 
as well as tree amplitudes. If a precise evaluation is needed (below 10% uncertainty), 
then the next-to-next-to-leading order terms, involving two-loop amplitudes, may also 
be required. 

In principle, Feynman rules |@] are all we need to evaluate the tree and loop ampli¬ 
tudes. In practice, however, although Feynman rules are very general, applying to any 
local quantum field theory, by the same token they are not optimized for the problems 
at hand. More efficient methods are available, which make use of the extra symmetries 
(some hidden) of QCD. 


TRANSFORMING TO TWISTOR SPACE 


An easy way to see that there should be more efficient methods out there is to notice 
that many QCD amplitudes are much simpler than expected. For example, the tree-level 
amplitudes for the scattering of n gluons turn out to all vanish, if the helicities of the 
gluons (considered as outgoing particles) are either a) all the same, or b) all the same, 
except for one of opposite helicity. Using parity, we can take the bulk of the gluons to 
have positive helicity, and write this vanishing relation as 

Ar(l=",2+,3+,...,n+)=0. (1) 


This vanishing is somewhat mysterious from the point of view of Feynman diagrams. On 
the other hand, it can be demonstrated simply using Ward identities arising from a secret 
supersymmetry that tree-level QCD amplitudes possess 0]- This symmetry allows two 
of the gluons to be replaced by their superpartners, gluinos, which can be taken to be 
massless here. Helicity conservation for the gluinos then implies the vanishing of the 
amplitudes. 

The first sequence of nonvanishing tree amplitudes has two gluons with negative 
helicity, labelled by j and /, say, and the rest of positive helicity. This sequence of 
maximally helicity-violating (MHV) amplitudes has an exceedingly simple form i,S], 
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in terms of spinor products (z j) we shall define shortly. Equation @ is the expression 
for, not the full amplitude, but rather a piece of it where the n gluons have a definite 
cyclic ordering. The full amplitude can be built out of permutations of such partial 
amplitudes, as reviewed for example in ref. |0]. Some of the structure of eq. @ follows 
from supersymmetry, but not all. 

To see much more of the structure, Witten JIl] transformed the amplitudes Q from 
the traditional momentum-space variables, into a twistor space invented by Penrose jH]. 
The twistor transform is a kind of Fourier transform. There are many examples where 
transforming a problem into the right variables can expose its simplicity. For example, 
if we measure the time dependence of the electric field E{t) associated with the light 



emerging from some glowing sample of gas, we find a fairly unenlightening waveform. 
However, if we use a speetrometer to measure instead the frequeney (energy) speetrum 
of the light, that is, E{g)) = / dte’^^E{t), we find spectral lines, which are clues toward 
decoding the structure of the emitting gas. In in an analogous way, the twistor trans¬ 
form exposes certain lines on which QCD amplitudes are localized or supported, thus 
revealing more of their structure, and giving rise to new, more efficient ways to compute 
them. 

Before describing the twistor transform, however, we should discuss the spinor vari¬ 
ables used in eq. Q, because they are well-suited for describing scattering amplitudes 
for massless particles with spin, and are the starting point for the twistor transform. Let 
i=l,2,...,n label the particles being scattered. Usually, the four-momentum vectors , 
which transform under the spin-1 representation of the Lorentz group, are used as the ar¬ 
guments of the amplitude, A = A{ki). The relativistic invariants constructed out of these 
vectors are the Lorentz inner products, or invariant masses, Sij = 2ki ■ kj = (ki + kj)^, 
which are equivalent in the massless case, kf = 0. However, for massless particles with 
spin, it is better to “take the square root” and use, instead of kf, objects transforming 
as the spin-1/2 representation of the Lorentz group, namely the massless Dirac spinors 
associated with momentum k,-, M±(k,), where the ± sign labels the helicity. A shorthand 
notation for the two-component (Weyl) versions of these spinors is, 

{Xi)a=[u+{ki)]^, {li)a=[u-{ki)]^. (3) 

We can always reconstruct the momenta from the spinors, using the positive-energy 
projector for massless spinors, u{k)u{k) = k, or in two-component notation, 

kf {(y^)aa = i^i)aa — ■ (4) 

Equation © shows that a massless momentum vector, written as a bi-spinor, is the 
product of a left-handed spinor with a right-handed one. 

Instead of Lorentz inner products of momenta, Sij = 2k, ■ kj, we use spinor products, 
defined by 

{jl) =fi_(k;)M+(k/), [;■/] = = u+{kj)u-{ki) , 

(5) 

where and are antisymmetric tensors for SU (2). These products satisfy 

(* j) [j = ^Tr[k,-k;] = 2ki ■kj = sij. (6) 

So they are just the square roots of the Lorentz inner products, up to a phase (/), 

{jl) = , [j 1] = . (7) 

The utility of spinor variables for QCD amplitudes was recognized already in the 
1980s 0,0]. They precisely capture the “square-root-plus-phase” behavior of gauge 
theory amplitudes as the momenta of two of the particles, i and j, become collinear. 
This behavior arises because the sum of the helicities of the final-state particles is never 




FIGURE 2. Contrasting collinear behavior of amplitudes in (a) massless scalar (j)^ theory, where there 
is no angular-momentum mismatch, and (b) massless gauge theory, for example a gluon splitting into two 
gluons, where there is always a mismatch. 


equal to the helieity of the almost-on-shell intermediate particle, as illustrated in fig.Etb) 
for the case of a gluon splitting into two gluons, for which ±1 ± 1 7 ^ ±1. This mismatch 
in angular momentum along the collinear direction lessens the singularity, from 1 /Sij 
(the behavior of the scalar theory shown in fig. Eta)) to 1 /h also introduces a 
phase depending on the azimuthal angle, which is conjugate to the angular momentum. 
Equation © shows that both characteristics are captured by putting a spinor product 
in the denominator of the amplitude, explaining why the spinor products are natural 
variables to use. In other words, we should write A = A,) instead of A = A{ki). 

Now we can describe the twistor transform i,[l . It is a “half” Fourier transform, in 
which the right-handed spinors A, are left untouched, but each left-handed spinor A,- is 
exchanged for its Fourier conjugate variable /i;, defined by 
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(These relations are completely analogous to the standard Fourier relation between mo¬ 
mentum and position, x = id/dp, p = —id/dx.) Since the spinors and their conjugates 
each have two components, twistor space has four coordinates (for each external par¬ 
ticle), (Ai,A 2 ,/i\/i^). However, because amplitudes are only defined up to a phase 
associated with external states, two points in twistor space are equivalent if the four 
coordinates differ by a constant multiple ^ (the complexification of the phase). 


(Ai,A 2 ,At\M^) = (^Ai,<^A 2 ,<^At\<^M^)- (9) 


So in fact (projective) twistor space is three-dimensional. 

What do amplitudes look like in this space? We can compute them by Fourier trans¬ 
forming just as we would to take a wave-function from position-space to momentum- 
space 111], 

A(A,-,A,) ^ A(A,-,M/)= /nA-e'>'%(A,-,A,0. (10) 

1=1 

The simplest cases to consider are the MHV amplitudes 0 , which contain only an¬ 
gle brackets ((/ j)), and so depend almost exclusively on the right-handed spinors A/, 
AMHV(a.^^.) Their only dependence on the left-handed spinors is through 

the usual momentum-conserving 5-function (which was implicit in eq. ©). This factor 







(a) MHV 


(b)NMHV 


(c) NNMHV 


FIGURE 3. Tree amplitudes for n gluons are supported on networks of intersecting lines in twistor 
space. The number of lines is one fewer than the number of negative-helicity gluons. 


can be written, using the identity 

5'^{k) = f d^xexp[ik-x] 


and eq. ©, as 


5"^ ( j = / d'^x exp ) a (^/) 
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Then the transformed amplitudes are 

{Xi,Hi) = /]^JA/exp[ijU/A,] /j'^j(;Ajf^^(A/)exp[/xA,A,] 

1=1 

= f d'^x(Xi) f Y[dh&w[Kl^i + xXi)Xi\ 

•' i=l 

= A^^'^{Xi) fd'^xYlSilH+x^i)- (13) 

i=l 

The product of all the linear 5-function constraints simply means that the amplitude is 
supported on a line in twistor space, as shown in fig.|3ta). 

More complicated amplitudes can also be inspected. The first nonvanishing, non- 
MHV n-gluon amplitudes are the six-gluon amplitudes with three positive and three 
negative helicities, first computed, from 220 Feynman diagrams, in 1988 |@]- The sim¬ 
plest case, where the three positive helicities are adjacent, is given by. 


Af‘=(l+,2+,3+,4^,5',6^) = i 
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where Sabc = (ka + kb + kc)^ and (a^| (b + c) \d^) = u-(ka)(kh + kc)u^ (k^Y 

The seven-gluon amplitudes were also eomputed around this time iflffl . using off- 
shell recursive methods to avoid dealing directly with the 2,485 Feynman diagrams. 
The explicit results in this case fill several pages. Computing the twistor transform via 
eq. (fTOt is rather difficult. However, suppose one has a guess for how the amplitudes are 
supported in twistor space, for example that they are localized on some curve described 
by a polynomial equation C(Z,) = 0, where Z = (Ai, A 2 , jU \ . Then it is relatively easy 
to check such a guess back in spinor-space, where C(Z,) becomes a differential operator, 
since /i/ = id jdXi. Applying C(Z,) to A(A/, A/), if the result vanishes identically then the 
amplitude is supported on the curve; that is, either C = 0 or else A = 0. 

This method was used last year by Cachazo, Svrcek and Witten Emmi to build up 
evidence for the picture illustrated in fig.|3l Scattering amplitudes for n gluons, of which 
have negative helicity, are localized in twistor space on networks of intersecting 
lines, where the number of intersecting lines is n- — 1. The MHV case, n- = 2, was 
discussed above. The next-to-MHV (NMHV) amplitudes with = 3, for example the 
six-gluon example in eq. d, are sums of terms, each of which is supported on a pair of 
intersecting lines, as shown in fig. |3tb). The partitioning of points among the lines can 
vary from term to term. Three intersecting lines are needed to describe the next-to-next- 
to-MHV (NNMHV) amplitudes with = 4 (fig. He)), and so on. 


MHV RULES 

While the twistor structure shown in fig. Elis extremely appealing, it does not directly 
yield the numerical values of the amplitudes. However, Cachazo, Svrcek and Witten ImAIi 
also wrote down a set of diagrammatic “MHV” rules, which can be used in place of 
Feynman rules to compute the amplitudes, and which make the twistor structure in 
fig. El manifest. Each MHV diagram generates a term in the amplitude which has one 
of the possible twistor structures, taking into account the possible partitionings of points 
among the (n_ — 1) lines. For example, the MHV diagram in fig. IH for an amplitude 
with = 4, corresponds to the twistor structure in fig. Etc). The helicities of internal, 
as well as external, gluons are labeled by ± in the diagram. Each vertex must have 
exactly two negative-helicity gluons attached to it, but it can have an arbitrary number 
of positive-helicity gluons, just like the MHV amplitude Q. In fact, the rule for this 
MHV vertex (the complex number associated with it) is given by eq. Q, with a simple 
prescription for continuing intermediate legs off shell. The rule for an internal line is 
a factor of ijp^, much like a scalar propagator. For processes with a large number of 
gluons, there are considerably fewer MHV diagrams than Feynman diagrams, because 
many Feynman subdiagrams get lumped into single MHV vertices. Also, the algebra 
required to evaluate each diagram is considerably simpler than for the typical Feynman 
diagram, because there is no tangle of Lorentz indices to follow. 

Because the MHV rules are so efficient, they were quickly generalized to a more 
general set of processes of interest in the context of LHC signals and backgrounds: tree- 
level QCD amplitudes containing massless external fermions as well as gluons Q; 
those with a Higgs boson, which couples to gluons via HTr{G^vG^'^) in the large 
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FIGURE 4. Example of an MHV diagram, corresponding to fig. 3(c). 


nit limit 111511 : and amplitudes includin g on e or more electroweak vector bosons in 
addition to massless quarks and gluons lll^ . A set of scalar-type rules for QCD with 
massive quarks {e.g. the top quark) was also produced, starting directly from the QCD 
Lagrangian id. 


TWISTOR STRUCTURE AT ONE LOOP 


In parallel with the extension of tree-level MHV rules to different processes, the twistor 
structure of one-loop amplitudes began to be investigated d. For multi-particle pro¬ 
cesses, one-loop amplitudes are much more intricate than tree amplitudes. Their twistor 
structure is also complicated by a “holomorphic anomaly” ifisll . in which derivatives 
from Hi = id/dXj act near singular regions of the loop integration. For these reasons, it 
has proven simpler to proceed by first representing amplitudes as linear combinations of 
various types of basic one-loop integrals — boxes, triangles, bubbles, etc. — and then 
examining the twistor structure of the coefficients of these integrals. 

The simplest situation to consider is a “toy model” for perturbative QCD, namely 
its maximally supersymmetric cousin, ^ = 4 super-Yang-Mills theory. In this theory, 
the coefficients of the triangle and bubble integrals all vanish, reducing the problem 
to that of determining the coefficients of box integrals d. These coefficients can be 
found quite readily lIlM 1^ d by inspecting either standard two-particle unitarity 
cuts or (more efficiently) generalized cuts ||25l dl where four propagators 

are held open UM- 

The resulting twistor structure lEfll |d Hzll is illustrated in fig. |5l In the MHV case 
shown in fig. |5ta), the only nonvanishing box coefficients are those where two of the 
external momenta for the scattering amplitude, and S 2 , are also momenta for the box 
integral; the remaining external momenta are partitioned into two diagonally opposite 
clusters, A and B. This integral is referred to as a two-mass box, because the clustered 
momenta Ka = h and Kg = 'Liesh are massive, The coefficient of the 

two-mass box lll 9ll is just the MHV tree amplitude (l2l), which is localized on a single 
line in twistor space (see fig. Ha)). In fig. |5ta), the single line has been redrawn as 
a pair of lines intersecting in two points, si and S 2 , to make its appearance consistent 
with an “MHV rules” approach to one-loop amplitudes ll^ . and with the pattern found 
for more negative-helicity gluons. (Just as in Euclidean space, a pair of straight lines 
intersecting in two points in twistor space is the same as a single line.) In the NMHV 
case shown in fig.|5tb), the simplest nonvanishing box coefficients are generically those 
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(c) NNMHV 


FIGURE 5. Twistor structure of box integral coefficients for one-loop amplitudes in jY = 4 supersym¬ 
metric Yang-Mills theory. For each series of amplitudes, (a) MHV (n_ = 2), (b) NMHV (n_ = 3), and (c) 
NNMHV (n_ = 4), the type of box integral having the simplest nonvanishing coefficient is depicted at the 
top, and the localization of those coefficients in twistor space is shown at the bottom. 


of the three-mass box integral, for whieh three of the legs, A, B, and C, represent 
elusters of momenta from the seattering amplitude, and only one, s, is an individual 
seattering momentum. These coefficients have a planar twistor structure, consisting of 
three intersecting lines, and the leg s sits at one of the intersections For 

the NNMHV case in fig. |5tc), the four-mass box coefficients have the nonplanar ring 
structure shown ll^ . In general, as in the tree case, fig. 0 one-loop box coefficients 
are supported on networks of lines, but the lines are joined into rings to match the 
loop topology. Similar structures have been found for coefficients of integrals in gauge 
theories with ^ < 4 supersymmetries |j^. 


WHAT IS A TWISTOR STRING? 

I have been remiss in titling this talk “Twistor String Theory and QCD,” without say¬ 
ing anything yet about what twistor string theory is, or how it is related to the more 
phenomenological developments just outlined. In fact, I won’t describe twistor string 
theory at any length, but I would like to briefly contrast it with ordinary string theory, 
from the perspective of methods for computing gauge theory amplitudes. 

An ordinary string is an extended object which moves in space-time. Different phys¬ 
ical vibrations of the string are associated with different particle states. The higher the 
harmonic, the more massive the particle; indeed, there is an infinite tower of ultra-heavy 
particles, as well as a set of massless ones. One of the massless particles is always the 
graviton. Because the one energy scale in gravity is the Planck mass, Mp Ri 10^^ GeV, 
this sets the scale for the ultra-heavy masses (unless certain extra dimensions happen to 
have a large size), as shown in fig. |3a). Also, the massless spectrum can be relatively 
complicated — several gauge groups, matter fields transforming in various ways, and so 
on. 
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FIGURE 6. (a) Typical spectrum of particles in ordinary string theory, (b) Typical spectrum in (topo¬ 

logical) twistor string theory. 


In the early 1990s, this type of string theory was adapted by Bern and Kosower 
into a tool to eompute one-loop QCD amplitudes a. It worked pretty well. The 
first eomputations of the helieity amplitudes for four-gluon seattering in QCD used 
this teehnique lUO]- (Th e unpolarized eross seetions were eomputed earlier via more 
traditional methods Bill . ) Th e five-gluon helieity amplitudes were also eomputed via 
the string-based method ll^ . On the other hand, this string theory was not optimized 
for QCD ealeulations. It did not possess all the symmetries of QCD. For example, QCD 
is elassieally eonformally invariant (independent of energy seale). But traditional string 
theory depends on the seale Mp, as refieeted in the mass speetrum in fig.l^^a). Related to 
this faet, a fair amount of analysis was required to deeouple the unwanted massive states 
from the loop amplitudes (as well as the massless states not eorresponding to QCD). 
The analysis would have had to be redone to ineorporate external quarks, for example. 
In the meantime it was found that “abstraeting the lessons” from string theory was often 
the most effieient way to proeeed. The effieieney of the string-based rules for amplitudes 
eould be attributed to a baekground-field quantization of gauge theory ll^ and a seeond- 
order formulation for fermions, for instanee a . These lessons eould then be applied 
to amplitudes with external quarks, without having to develop the full string-theoretie 
maehinery li^ . 

In eontrast, the twistor string theory invented by Witten lUl] is a topologieal one, 
and the string moves in twistor-spaee, not the usual spaee-time. “Topologieal” means 
that the energy of the string only depends on topologieal information, so that very 
few of its degrees of freedom are dynamieal. As a result, it does not have a tower 
of massive states, only massless ones, as shown in fig. |6tb). Twistor string theory is 
eonformally invariant, like elassieal QCD. It makes mueh more manifest the symmetries 
of elassieal QCD, whieh inelude not only eonformal invarianee, and the seeret (^ = 4) 
supersymmetry mentioned in seetion„ but a full supereonformal group eontaining them. 
So, from the point of view of ealeulating QCD amplitudes, twistor string theory seems 
almost designed to do the job. 

On the other hand, twistor string theory is still not preeisely QCD. It possesses all 
the jV = 4 superpartners of the gluons, instead of quarks. It also eontains gravitons, 
but not those of Einstein’s theory of gravity; instead they belong to a non-unitary theory, 
eonformal supergravity. Both of these properties are not really an issue for eomputations 







of tree-level amplitudes, but they ean play havoe with a loop-level deseription. In faet, 
there is no satisfaetory one-loop formulation of twistor string theory at present. Onee 
again, however, from a computational point of view, abstracting the lessons is often the 
best. 

Even at tree level, such abstraction can be beneficial. Althou^ the MHV n-gluon 
amplitudes Q could be evaluated directly from the twistor string |Il[|, and the six-gluon 
non-MHV amplitudes, such as eq. (fT4l) . were also produced in this way ll3^ . the MHV 
rules d have provided a much more efficient method for generic tree amplitudes. 
They originated at least in part from abstracting the twistor structure which was found 
by studying existing QCD amplitudes. (One could also say, however, that the MHV rules 
follow from a different, “disconnected”, prescription for evaluating the relevant twistor- 
string contributions.) 


ON-SHELL RECURSIVE METHODS 


Another process of abstraction and streamlining led, at the beginnin g o f this year, to the 
on-shell recursion relations of Britto, Cachazo, Feng and Witten 11^ 1^ . These relations 
are even more efficient, and lead to more compact formulas, than the MHV rules. Also, 
they can be proven in a very simple way, using only Cauchy’s theorem and factorization 
properties. So it is very easy to extend these relations to more general processes, and 
also to apply the same kinds of techniques to the computation of one-loop amplitudes in 
QCD. 

The path to the on-shell recursion relations was somewhat roundabout, proceeding 
through the one-loop amplitudes in JV = 4 super-Yang-Mills theory, whose box co¬ 
efficients were sketched in section „ These amplitudes have infrared divergences, repre¬ 
sented in dimensional regularization as poles in e = (4 — D)/2. The residues of the poles 
have to be proportional to the corresponding tree amplitude. This req^uirement gave new 
formulas for tree amplitudes, in terms of sums of box coefficients 112 (1.1^1,3^ . which 
were more compact than previously-known expressions. Using generalized unitarity, 
these formulas could be reinterpreted as quadratic recursion relations jj^ 

The basic on-shell recursion relation for tree amplitudes reads ll^ l38ll. 


Ar(l,2,...,n) 
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(15) 


It is depicted diagrammatically in fig. |71 The amplitude is represented as a sum of 
products of lower-point amplitudes, evaluated on shell, but for complex, shifted values 
of the momenta (see below). The helicity labels of the n external gluons have been 
omitted, but they are the same on the left- and right-hand sides of eq. (fTSb . For the 
relation to be valid, the helicities of gluons n and 1 can be (fi„,fii) = (—1,1), (1,1), 
or (-1,-1), but not (1,-1). There are two sums. The first is over the helicity h of an 
intermediate gluon propagating (downward) between the two amplitudes. The second 
sum is over an integer k, which labels the different ways the set {l,2 ,...,n} can be 
partitioned into two cyclicly-consecutive sets, each containing at least 3 elements, where 



n-1 



FIGURE 7. Diagrammatic representation of an on-shell recursion relation for tree amplitudes. 


labels 1 and n belong to different sets. A hat on top of a momentum label denotes that 
the eorresponding momentum is not that of the original n-point amplitude, but is shifted 
to a different value. 

To deseribe the shifted momenta, first note that, from eq. is a 

t \JL tX t L 

massless four-veetor beeause of the antisymmetry of the spinor produets. 


tj = = - (;,) [i,] = 0. 


ai $ 


(16) 


It will eontinue to be massless even if one of the two spinors is shifted so that it is no 
longer the eomplex eonjugate of the other spinor, for example 


(<7'p)ao! — (^i)aa — (^!)a(^z)o! ? 
where A,- is shifted away from A, . 

The momentum shift in the term in eq. (fT51) ean now be deseribed as. 


(17) 


Ai 

A„ 


Ai = Ai +Zk^m 

A,,, 


Ai Ai, 


A,j — A^' 


■ZfcAi, 


where 


Zk 




{n-\^l,kn~) 


(18) 


(19) 


This shift keeps = (Ai + Zk^n)^i and = A„(A„ —Zy^Ai) massless, as diseussed 
above. It preserves overall momentum eonservation, beeause ki+kn = AiAi + A„A„ = 
k\+kn. And the intermediate gluon momentum, defined by iti ^ = iti + z^A„Ai, is also 
massless (on shell), beeause 

Kl, = {li^k + ZkW^ = Kl, + Zk{n-\li,k\l-) = 0. ( 20 ) 

The derivation of eq. ([131 is very simple ll^ . The momentum shift (fTRll is eonsidered 
for an arbitrary eomplex number z, instead of the diserete values Zk in eq. dlSll. This 




shift defines an analytie funetion of z, It has poles in z whenever a eollection 

of the shifted momenta, eorresponding to an intermediate state, can go on shell. For 
every allowed partition of {1,2,..., n} into {1,2,. .., U + 1,..., n — 1, n}, there is a 
unique value of z that accomplishes this, Zk, because ^(z*) = 0 according to eq. (l 2 (H . 
The desired amplitude is the value of A5f®®(z) at z = 0. Provided that > 0 as 

z —oo, this value at z = 0 is determined by Cauchy’s theorem in terms of the residues 
of at z = Zk- Using general factorization properties of tree amplitudes, the 

residue evaluates to the product found in the term in eq. (fTSt . The vanishing of 
as z ^ oo can be established directly from Feynman diagrams for (fi„, hi) = ( —1,1) 11^ . 
For the other two valid cases it can be shown using the “MHV rules” lll^ . or by a 
recursive argument S]. 

No knowledge of twistor space is needed to implement eq. O. Its derivation is 
heuristically related to twistor space, however, in that spinors, not vectors, play the 
fundamental role. 

Off-shell recursive approaches to summing Feynman diagrams have a long his¬ 
tory lEHliil]. In the off-shell case, however, the auxiliary lower-point quantities are 
gauge-dependent. In on-shell recursion relations, in contrast, they are precisely the de¬ 
sired physical, gauge-invariant, on-shell scattering amplitudes, just with fewer partons. 
In short, trees are recycled into trees. 


A simple application: Ag®®(l+,2+,3+,4 ,5 ,6 ) 

Let us now work through a simple application of eq. (fTSll ll^ . The first non- 
MHV n-gluon amplitudes (taking into account parity) are those with six gluons, 
three of positive helicity and three negative. There are three cyclicly-inequivalent 
helicity configurations: Ag®®(l+,2+,3+,4^,5^,6^), Ag®®(l+,2+,3^,4+,5^,6^), 
and 2^,3+,4^,5+, 6 ^). The last of these ai^litudes is related to the first 

two a “dual Ward identity” (group theory relation) Bj]. Here we apply eq. (fTSl) to 
Ag®®(l + ,2+,3+,4^,5^,6^). Instead of 220 Feynman diagrams (including all color- 
orderings), there are just three potential on-shell recursive diagrams, shown in fig. [5] 
Diagrams of the form of fig. [Hta) and fig. [Etc), but with a reversed helicity assign¬ 
ment to the intermediate gluon, vanish because A 3 ®®(-|-,-|-,-|-) = A 3 ®®(—,—,—) = 0 . 
Figure IHtb), and the corresponding diagram with a reversed intermediate helicity, both 
vanish using eq. dU)- Finally, diagram (c) is related to diagram (a) by the “flip” symmetry 
(l<-^6,2'«-»-5,3'«-»-4) (plus spinor conjugation). 

So, remarkably, there is only one independent diagram, fig. Eta). Its value is given 
by the product of two shifted MHV amplitudes, each parity-conjugated with respect to 
eq. (El), 


D(a) ^ Ar(i + ,2+,-42)^<‘^(^t2^3+,4-,5-,6-) 

^ -i [ 12 ]^ [^ 1,2 3 ]^ 

^12 [2^1,2] [^1,2 1] [34] [45] [56] [6^1,2] 







(a) (b) (c) 


FIGURE 8. On-shell recursive diagrams for (l + ,2+,3+,4 ,5 ,6 ). 


(6“|(1+2)|3- 
'(61)(12)[34][45]^612(2- 

To get the eontribution , we add the image of under the permutation (l<-^6,2'«-»- 
5,3 4), eombined with spinor eonjugation, () ^ [ ]. The full amplitude is 




(6+l)|5->- 


( 21 ) 


Ar(l+,2+,3+,4-,5-,6-) 


(6-|(l + 2)|3-)^ 

(61)(12) [34] [45].612(2-|(6 + 1)|5-) 
(4-|(5 + 6)|r)^ 

(23)(34) [56] [61]556i(2-|(6+1)|5-). 


( 22 ) 


Let’s eompare this representation of the amplitude with the previous one, eq. d, found 
using Feynman diagrams. The second expression is shorter. (Here the difference in 
length is minimal; it becomes more striking for seven gluons lIlflLl^ .l It also makes 
manifest the square-root collinear behavior in all channels. For example, in the collinear 
limit where becomes parallel to ^ 4 , eq. (1^ has the correct 1/ (34) and l/[34] 
behavior manifest; in eq. (O, cancellations between the three terms, each of which 
behaves like I/ 534 , are required to obtain the proper behavior. On the other hand, eq. (1^ 
contains a spurious singularity, because ( 2^|(6 + 1)|5^) vanishes when k(,-\-ki happens 
to be a linear combination of ^2 and k^ (use the massless Dirac equation to see this). 
The amplitude is perfectly finite in this region, but each term diverges. In numerically 
implementing eq. (l 2 ^ . one should take care in this region. 


ONE-LOOP AMPLITUDES 

Although the MHV and on-shell recursive rules are quite efficient for the analytical 
computation of many types of tree amplitudes, and shed a lot of light on their structure, 
in the end all one really wants are numerical values. Quite efficient numerical computer 
programs have already been developed over the years, based on off-shell recursive 
methods il, which can evaluate QCD tree amplitudes with of order 10 external partons 










in a reasonable amount of time. In contrast, the complete set of one-loop helicity 
amplitudes is not known for any pure QCD process with greater than five external legs. 
There are similar bottlenecks for processes in which a few electroweak vector bosons are 
produced in addition to multiple QCD partons. So it is of great interest to see whether 
new methods can be developed for one-loop QCD amplitudes. 

The method used to prove the tree-level on-shell recursion relations ll^ — shifting 
a pair of momenta by a complex amount, while keeping them on shell — is particularly 
promising in this regard, because it efficiently incorporates the known factorization of 
amplitudes onto collinear and multi-particle poles. Indeed, the same techniques can 
be adopted at one loop, in order to determine the rational (non-logarithmic) parts of 
amplitudes, once the parts containing branch cuts (logarithms, poly logarithms, etc.) 
have been determined by other means — for example, using unitarity, as mentioned in 
section „ Recently, all the one-loop n-gluon helicity amplitudes in QCD with up to two 
(adjacent) negative-helicity gluons, and an arbitrary number of positive helicity ones, 
have been produced (or reproduced) in this way 114^14^ . The amplitudes having n =0 
or 1 are quite special, because they vanish at tree-level (eq. ([T]))- They have no infrared or 
ultraviolet divergences, and there are no branch cuts at all. Also, they were known from 
previous work 114511^ . They have a purely recursive representation, whose construction 
involved a few assumptions, which could be cross-checked by comparing to the previous 
results a. 

The series of amplitudes 3+,4+,...,«+), with two adjacent negative 

helicities, have branch cuts as well as infrared and ultraviolet divergences. The branch 
cuts were determined a decade ago, using unitarity ll^ . The rational parts can now 
be constructed recursively a . In addition to a set of recursive diagrams, much like 
the tree-level formula (El), there are certain “overlap” diagrams, which perform book¬ 
keeping with respect to certain rational-function terms which naturally accompany the 
logarithmic terms. There are relatively few diagrams to evaluate. For the rational part of 
Ag~i°°P(l^,2^,3+,4+, 5+, 6+), for example, there are four nonvanishing recursive dia¬ 
grams and three nonvanishing overlap diagrams. The evaluation of each diagram is com¬ 
pletely algebraic; no loop integrations are required. In contrast, the number of one-loop 
6 -gluon Feynman diagrams in QCD is 10,680, each of which requires a loop integration. 


CONCLUSIONS 

Several new methods for computing gauge theory scattering amplitudes relevant for 
LHC physics have been developed over the last year or two, with a strong stimulation 
from twistor string theory. After some abstraction and streamlining, however, many of 
these methods actually bear a close resemblance to the bootstrap program developed 
in the 1960s. In a bootstrap, scattering amplitudes are reconstructed directly from their 
analytic properties, without the need for a Lagrangian ll^ l^. While this program has 
proven difficult, if not impossible, to carry out in full nonperturbative generality in a 
strongly-coupled four-dimensional field theory, in the context of perturbation theory 
much more information is available to assist it. The (factorization) poles of amplitudes 
are dictated by amplitudes with fewer legs, while the (unitarity) cuts are dictated by 


products of amplitudes with fewer loops. Tree-level on-shell reeursion relations, for 
example, are a very eonvenient way of systematieally ineorporating the faetorization 
data. The use of analyticity fell somewhat out of favor in the the 1970s, with the rise 
of a Lagrangian (QCD) for the strong interaetions. Ironieally, it now proves useful to 
resurreet analytieity, and a perturbative bootstrap, as a tool for eomputing eomplieated 
QCD amplitudes — for which a direct Lagrangian approach, that is, using Feynman 
rules, ean be very cumbersome. 

To date, the “praetieal” spinoffs from twistor-inspired methods have been primarily 
for tree amplitudes (whieh ean also be obtained by other, numerieal methods), and 
for loop amplitudes in supersymmetrie theories. But reeently, new one-loop helicity 
amplitudes in full QCD have begun to fall to these methods, suggesting that soon there 
will be direct phenomenological applications. In addition, the recent rapid progress in 
developing new eomputational approaehes along these lines is very likely to eontinue. 
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